High cumulants of current fluctuations 
out of equilibrium 
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We consider high order current cumulants in disordered systems out of equilibrium. They are 
interesting and reveal information which is not easily exposed by the traditional shot noise. Despite 
the fact that the dynamics of the electrons is classical, the standard kinetic theory of fluctuations 
needs to be modified to account for those cumulants. We perform a quantum-mechanical calculation 
using the Keldysh technique and analyze its relation to the quasi classical Boltzmann-Langevin 
scheme. We also consider the effect of inelastic scattering. Strong electron-phonon scattering renders 
the current fluctuations Gaussian, completely suppressing the n > 2 cumulants. Under strong 
electron-electron scattering the current fluctuations remain non-Gaussian. 



I. INTRODUCTION 

The fact that electric current exhibits time dependent 
fluctuations has been known since the early of 20 th cen- 
tury. Still it remains an active field of experimental 
and theoretical research 1 . Among other features, non- 
equilibrium shot noise can teach us about the rich many- 
body physics of the electrons, and may serve as a tool 
to determine the effective charge of the elementary car- 
riers. While the problem of non-interacting electrons is 
practically resolved by now, the issue of noise in systems 
of interacting electrons remains largely an open prob- 
lem. In this paper we focus on two main issues: (1) 
The role of both inelastic electron-electron and electron- 
phonon scattering for current correlations. (2) In the 
absence of interference effects one is tempted to employ 
the semi-classical kinetic theory of fluctuations 2 ' 3 , also 
known as the Boltzmann-Langevin scheme. This ap- 
proach has been originally developed to study pair corre- 
lation functions, hence is not naturally devised for higher 
order cumulants. A naive application of this approach 
turns out problematic. To study high order cumulants, 
we first need a reliable scheme, which is why we resort to 
a microscopic quantum mechanical approach. The com- 
parison between microscopic calculations and the appli- 
cation of the Boltzmann-Langevin scheme for the study 
of higher order cumulants is the second objective of this 
paper. 

The outline of this work is the following: Section I 
addresses a few introductory issues concerning current 
fluctuations in non-interacting systems. These brief re- 
minders are necessary for the following analysis. In Sec- 
tion I A we consider the two main sources of noise in non- 
interacting systems: thermal fluctuations in the contact 
reservoirs and the stochasticity of the elastic scattering 
process involved. For this purpose a simple stochastic 
model is studied. In Section IB we explain why the 
study of high order correlation functions is of interest. 
In Section I C we recall some elements of the kinetic 
theory of fluctuations, explaining the difficulty in apply- 



ing it to high cumulants. Section II is devoted to the 
analysis of the third order current cumulant. We con- 
sider two limiting cases: that of independent electrons 
(II A), and that of high electron-electron collision rate 
(II B). In Section III we compare the results of the mi- 
croscopic quantum mechanical analysis with those of the 
Boltzmann-Langevin scheme. In Section IV we discuss 
higher moments of current cumulants the full counting 
statistics, and comment briefly on the effect of electron- 
phonon scattering. 



A. Noise In Non-Interacting Systems: Probabilistic 
Scattering And Thermal Fluctuations 



Before discussing the problem of interacting electrons 
we would like to recall some important features of fluc- 
tuations in a non-interacting electron gas (for a review 
see 1 ) . One of the issues addressed in the present work is 
whether what has become to be known as 11 quantum noise 
" can be properly discussed within the (semi)-classical 
framework of the kinetic theory of fluctuations. Surely, 
non-equilibrium shot noise depends on the discreteness 
of the elementary charge carriers (and this charge is 
quantized). Also, at low temperatures one is required 
to employ the Fermi-Dirac statistics governing the occu- 
pation of the reservoir states. In addition, the channel 
transmission and reflection probabilities are governed by 
quantum mechanics. But other than that, interference ef- 
fects appear to play a very minor role in the formation of 
current fluctuations. Hence the appeal of a semi-classical 
approach. To understand the origin of the current fluc- 
tuations we consider a simple stochastic model, following 
earlier works 4-6 . This model is a caricature of the physics 
underlying low frequency transport through a quantum 
point contact (QPC) with a single channel having a trans- 
mission probability T. The QPC is connected to two 
reservoirs with respective chemical potentials fj,i and (j,r 
(this model is readily generalized to the multi-channel case) 
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FIG. 1. A simplified view of transport through a single channel quantum point contact. For definitions consult the text. 



For our needs electrons occupying (in principle broad- 
ened) single particle levels (depicted in a Fig.l) can 
be perceived as classical particles attempting to pass 
through the barrier (with transmission probability T). 
In principle T may be a function of the level's energy 
(T — > %). For the mean level spacing Ae the time 
At = h/Ae can be interpreted as a time interval between 
two consecutive attempts. Since the electron charge is 
discrete (and is well-localized on either the left or the 
right side of the barrier) the results of a transmission 
attempts have a binomial statistics. 

Consider the fluctuation in the net charge transmitted 
through the QPC over the time interval t 3> At. The 
number of attempts N — t/ At occurring within this in- 
terval is large. Since the attempts are discrete events, 
they can be enumerated. To describe the result of the 
n th attempt of an electron originating from level i on the 
left (right) we define the quantity Pi(n) and qi(n) (cf. 
Fig.l). For an attempt ending up in the transmission 
of an electron from left to right (from right to left) the 
quantity pi (n) (qi(n)) is assigned the value 1. Otherwise 
Pi(n) (qi(n)) are zero. This implies that {pi(n)}'s and 
{<li( n )Y s have each binomial statistics with the expecta- 
tion value {pi(n)} = {qi(n)} — %. We also assume that 
neither different attempts from the same level, nor at- 
tempts of the electrons coming from different levels are 
correlated. In addition to fluctuations in the transmis- 
sion process, there are also fluctuations of the occupation 
numbers (bi) of a single particle level i. These repre- 
sent thermal fluctuations in the reservoirs. For fermions 
they possess a binomial distribution 7 . We end up with 



P(x) = (1 - (x))S(x) + (x)S(l - x) 
for x =p i (n),q i (n),a i ,b i , 



and the correlations 



(Pi{n)pj(n')) = 8(i,j)5(n,n')%; 
(qi(n)qj(n')) = 5(i,j)8(n,n')%. 



(1) 



(2) 



i = 3 
i ^ 3- 



(bibj) 



o. 



l = 3 
i ^ 3- 



(3) 



Here nf = (ai) and nf = (bi). The typical time scale 
over which the occupation number of a level fluctuates is 
usually dictated by the interaction among the electrons 
or with external agents. We assume that at the leads the 
time fluctuations of the occupation numbers (ai and bi) 
over the interval t are negligible. 

Motivated by this picture, we consider fluctuations in 
the net number Qt of electrons transmitted through the 
constriction within the time interval t. Employing the 
Pauli exclusion principle one can write 



N 



Qt = Y2,\Pi( n ) a i{ 1 - b i) - Qi(n)bi(l - di)] 



• (4) 



According to eqs.(2,3 and 4) the expectation value of 
the transmitted charge is 



<Q t ->=TJV5>f(l-nf)-n?(l 



(5) 



Evidently it is related to the value of the d.c. current 
through (Qt) — il/e. We next discuss the higher or- 
der moments of current fluctuations. As was shown in 
Refs. 8 (see also 9 ) the experimentally measured high or- 
der current cumulants should be defined in quite a sub- 
tle manner. To make contact between the measured ob- 
servables and theoretically calculated quantities, we re- 
call that within the Keldysh formalism 10 the time axes 
is folded. For any given moment of time there are two 
current operators, one for the upper and another for the 
lower branch of the contour. The product of the symmet- 
ric combination of these two operators, I 2 , time ordered 
along the Keldysh contour (T c ) and averaged with respect 
to the density matrix yields the proper current correla- 
tion function. In a stationary situation the pair current 
correlation function 11 



Also 



S 2 (t-t') = (T c 7 2 (t)7 2 (t / )) 



(6) 
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depends only on the difference between the time indices. 
Similarly we define the third order current correlation 
function 

S 3 (t-t';t' -t") = {T c h{t)h{t')I{t")) . (7) 
In Fourier space it can be represented as 

S 2 (wi) - fd(t- t')e^ l(t -°S 2 (t - t') , (8) 

J c 

s^.^Jd^-tyit'-ty- 1 ^ 1 ^^'- 1 ''^ -t',t' - 1"). 

Usually it is more interesting to consider only the cumu- 
lants, i.e. irreducible parts of the correlation functions. 
For the pair current correlation function we have 

S 2 (lu = 0) = S 2 -I 2 , (9) 

while the third order current cumulant is given by 

5 3 (wi - 0, uj 2 = 0) = S 3 - 3/S 2 + 2/ 3 , (10) 

where S 2 ,S3 are taken at zero frequencies. Microscopic 
calculations of the current cumulants will be performed 
in the Sections II and IV below. Here we evaluate the 
cumulants within the stochastic model described above. 
The results agree with the low frequency current cumu- 
lants of the non-interacting electrons in the QPC 4 . The 
temperature dependence of those cumulants is qualita- 
tively similar to the dependence of the current cumulants 
in a disordered junction, analyzed below. 

To find the variance of the stochastic variable Q we 
employ eqs.(2, 3 and 4) 

(Q t 2 > - (Qt) 2 = T(l - T) £>i(i) - n 2 R (i)} 2 + 

i 

+ n R (i) - n£(i) - n 2 R (i)\. (11) 

i 

The first term on the r.h.s. of eq.(ll) vanishes at ther- 
mal equilibrium. We associate this part with the shot 
noise of the electrons. The second term vanishes at zero 
temperature, and we associate it with the thermal noise. 
As expected the "shot noise" part vanishes in the limit of 
perfect conductor (T — > 1). Going along the same pro- 
cedure for the third order cumulant one obtains 5 in the 
limit of low temperature (eV » T) 

((Qt- (Qt)f) = T(l - T)(l - 2T)<Q f > (12) 

and in the limit of high temperature (eV <C T) 

((Qi-(Qi)f)=T(l-T)(Q i ). (13) 

Note that with increasing the temperature the third or- 
der cumulant of the transmitted charge approaches a con- 
stant. This is a robust feature of all odd order cumulants, 
which can be understood from quite general arguments. 
Indeed, since the current operator changes sign under 



time reversal transformation, any even-order correlation 
function of the current fluctuations (e.g. <S 2 ) taken at 
zero frequency is invariant under this operation. Assum- 
ing that current correlators are functions of the average 
current, /, it follows that even-order correlation func- 
tions depend only on the absolute value of the electric 
current (and are independent of the direction of the cur- 
rent). In the Ohmic regime this means that even-order 
current correlation functions (at zero frequency) are even 
functions of the applied voltage. Evidently, this general 
observation agrees with the result eq. (11). 

By contrast, odd-order current correlation functions 
change their sign under time reversal transformation. In 
other words, such correlation functions depend on the 
direction of the current, and not only on its absolute 
value. Therefore in the Ohmic regime, odd-order corre- 
lation functions of current are odd-order functions of the 
applied voltage. This condition automatically guarantees 
that odd-order correlation functions vanish at thermal 
equilibrium. 

One can show that by considering high order moments 
of the stochastic model one reproduces the correct re- 
sults for any cumulants of a current noise in a QPC. 
Of course, the solution of a toy model can not replace 
the real microscopic calculations. But the fact that the 
results of the latter and the stochastic model agree sug- 
gests that the underlying physics is rather simple (and it 
is basically captured by such a simple model). It is the 
combination of thermal fluctuations (in the occupation 
of the single electron states) and random transmission of 
particles through the barrier that gives rise to current 
fluctuations. These two sources of stochasticity remain 
there when the electron can no longer be considered non- 
interacting. In that case, though, one cannot consider 
fluctuations at different energy levels to be independent. 



B. Why Are High Order Cumulants Interesting? 

Low frequency current fluctuations give rise to a large 
(in general infinite) number of irreducible correlation 
functions (cumulants). The pair current correlation func- 
tion provides us with only partial information about cur- 
rent fluctuations. To obtain the complete picture one 
should consider high order cumulants as well. As we have 
explained in Section I A the symmetry-dictated proper- 
ties of odd and even correlation function are very dif- 
ferent from each other; in particular odd order current 
cumulants are not masked by thermal fluctuations. For 
this reason they can be used for probing non-equilibrium 
properties at relatively high temperatures. Shot noise has 
been used to detect an effective quasi-particle charges in 
the FQHE regime 12,13 . Potentially, odd moments can 
be used to measure the effective quasiparticle charge in 
other strongly-correlated systems 6 . This may be needed 
for systems undergoing a transition controlled by temper- 
ature (for example normal-super-conductor) and having 
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different fundamental excitations at different tempera- 
ture regimes. One may hope therefore, that better un- 
derstanding of current correlations will teach us more 
about the many-body electron physics. At this moment 
this remains a challenge. Before trying to reach this goal, 
we need to understand the genuine properties of the high 
order cumulant in the relatively simple physical models. 
This is done next. 



C. Background and Issues to Be Discussed 

As was mentioned above it is quite appealing to try 
to discuss current noise in terms of the semi-classical ki- 
netic equation. To be more specific, let us consider a 
disordered metallic constriction. Its length L is much 
greater than the elastic mean free path The disorder 
inside the constriction is short-ranged, weak and uncor- 
rected. Under these conditions the electrons kinetics can 
be described by the Boltzmann equation: 

r/(p,r,t) = Col{f}. (14) 

Here Col{f} is the collision integral and 

t= (^+v-V r + eE-V p ) (15) 

is the Liouvillc operator of a particle moving in an ex- 
ternal electric field E. The pair correlation function of 
any macroscopic quantity may be found from the kinetic 
theory of fluctuations. Within this theory the distribu- 
tion function (/) consists of the coarse-grained (/) and 
fluctuating ( Sf) parts 

/(p,r,f) = /(p,r,i)+J/(p,r,t). (16) 

Since / is a macroscopic quantity (a quantity associated 
with a large number of particles) it must satisfy the On- 
sager's regression hypothesis (see Ref. 14 , Section 19). 

To sketch this hypothesis for the case of interest we 
consider a perturbation of the equilibrium distribution 
function 5f(p,r,t), small but substantially larger than 
the typical fluctuations of the distribution function. It 
follows then that the distribution function (with a high 
probability) will evolve toward the equilibrium state. Its 
relaxation dynamics is governed by the Boltzmann equa- 
tion and because the perturbation is small, the collision 
integral can be linearized. According to Onsager the cor- 
relation function of any macroscopic quantity, and in par- 
ticular (Sf(p, r, t)/(p',r',t')), is governed by the same 
equation as the one governing its relaxation (i.e. as equa- 
tion governing the quantity Sf(p, r, t)) 

(l(p, r, t) + J(p, r, *)) (5f(p, r, t)5f(p', r', t')) = 0, 

for t > t'. (17) 

Here X is a linearized collision integral. It was later sug- 
gested by Lax 15 (and can be proven within Kcldysh for- 
malism) that eq.(17) does hold for any stationary, not 



necessarily equilibrium, state. However, we need to re- 
call that the Onsager hypothesis was formulated only for 
a pair correlation function. There is no obvious way to 
apply this logic for higher cumulants. 

An alternative route of describing fluctuations within 
kinetic theory which is seemingly free of this difficulty 
had been proposed by Kogan and Shul'man 2 . Their 
picture is the following. The real space is divided into 
small volumes (as explained in a Ref. 14 ). The function 
/(p, r, t) represents the average number of particles in the 
state p of a unit volume element (cell) labeled by index 
(r) . The total number of the electron in every cell must 
be large. It was suggested in Ref. 2 that fluctuations of 
this number can be taken into account by adding a ran- 
dom (Langevin) source term to the Boltzmann equation. 
The resulting stochastic equation (including this additive 
noise) is called the Boltzmann-Langevin equation. 

(£(P, r, t) + Z(p, r, tj) 5 f(p, r, t) = 5J(p, r, t) . (18) 

The Langevin source 6 J(p, r, t) denotes a random num- 
ber of particles incoming into the given state in some 
interval (around time t). Since eq.(18) is a linear one, 
the statistics of the distribution function is determined 
by the random source term 8 J. To establish its proper- 
ties Kogan and Shul'man had used a rather simple phys- 
ical picture. To be consistent with the Boltzmann equa- 
tion they have assumed that interference effects are weak. 
The collision events are local in space and time. Since 
the typical number of electrons inside the cell is large one 
can ignore the correlation between the scattering of dif- 
ferent electrons. The electron scattering is a Poissonian 
process, with the number of scattered particles within 
any given cell (over a microscopic time interval) being 
large. While this picture yields a correct result for the 
pair correlation function, it substantially underestimates 
all high order correlators (starting from S3). 



II. THIRD ORDER CURRENT CUMULANT 

In this section we use microscopical calculations to 
evaluate the third order current cumulant for a quasi 
one-dimensional system of a length L with diffusive 
disorder 16 . We start with the coordinate-dependent cor- 
relation function 

S 3 (x,t;x\t';x",t") = (TJ 2 (x,t)i 2 (x\t')i 2 (x'\t")) . (19) 

Here x is a coordinate measured along a quasi one- 
dimensional wire (0 < x < L) of cross-section A; T c 
is the time ordering operator along the Keldysh contour. 
Next we perform the Fourier transform with respect to 
the time difference as in eq.(??). For small values of the 
frequencies uii,u>2, (small compared with the inverse dif- 
fusion time along the wire), the current fluctuations are 
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independent of the spatial coordinate. We next evalu- 
ate the expression, eq. (19), for a disordered junction, in 
the hope that the qualitative properties we are after are 
not strongly system dependent. In the present section 
we consider non-interacting electrons in the presence of 
a short-range, delta correlated and weak disorder poten- 
tial (e/r 3> fi, where r is the elastic mean free time and 
e/ is the Fermi energy). To calculate S3 we employ the 
tr-model formalism, recently put forward for dealing with 
non-equilibrium diffusive systems (for details see Ref. ). 
The disorder potential is (5-correlated: 



and 



(U dis (r) U dis (r')) 



5(t-t>), 



(20) 



where v is the density of states at the Fermi energy. 
The Hamiltonian we are concerned with is: 

H = H + H[ nt . 



(21) 



The motion of electrons in the disorder potential is de- 
scribed by: 



H = f 



Volume 



— (V - za) 2 + C/ di[ 
2m 



#(r). (22) 



Here ca/e is a vector potential. The Coulomb interaction 
among the electrons is described by 

H int = i J rfrrfr'*(r)*(r')Vb(r - r')*(r)*(r / ) , (23) 



where 



^o(r-r') = 



(24) 



A. Weak Inelastic Collisions 

Following the procedure outlined in Ref. 17 , we intro- 
duce a generating functional and average it over disor- 
der. Next we perform a Hubbard-Stratonovich trans- 
formation, integrating out fermionic degrees of freedom. 
Employing the diffusive approximation one obtains an ef- 
fective generating functional expressed as a path integral 
over a bosonic matrix field Q 

Z[a] = J VQexp(iS[Q,a\) . (25) 

Here the integration is performed over the manifold 

Q(x, t, h)Q{x, h,t')dti = S(t - t'), (26) 



the effective action is given by 

iS[Q, a] = -^Tr{D (VQ + i[a al a , Q]f + 4ieQ} , (27) 



7i 











i 


1 — 1 


, 72= [\ 


1 



(28) 



Tr represents summation over all spatio-temporal and 
Keldysh components. Here ai and &2 are the Keldysh 
rotated classical and quantum components of a. Here- 
after we focus our attention on ai, &2, the components 
in the direction along the wire. The third order current 
correlator may now be expressed as functional differen- 
tiation of the generating functional Z[a] with respect to 
a 2 



Ss^l—t 2 ,t2 — t 3 )- 



S 3 Z[a] 



8 Sa2{x 1 ,ti)Sa2(x2,t2)Sa 2 (x 3 ,t 3 ) 



■ (29) 



Performing this functional differentiation one obtains the 
following result 



S3^1 — *2j ^2 — £3) = 



16 



M(xi,ti)i (X2,x 3 ,t 2 ,t 3 ) + 



(xi,ti +->x 3l t 3 ) + (xi,ti ^x 2 ,t 2 ) + 
^^-M(x 1 ,ti)M(x2,t 2 )M(x 3 ,t 3 ) ). 



(30) 

Here we have defined 

I D (x,x',t,t') = Tr K <Q x ,t,tH2Qx',t',a2 ~ S t ,t'li \$x,x', (31) 



M(x,t) =Tv K ^jdt 1 {[Q x , tM ;V] Q xM tji^. (32) 

We employ the notations Q(x,t,t') = Q x ,t,t'\ Tr K is the 
trace taken with respect to the Keldysh indices; () de- 
notes a quantum-mechanical expectation value. The ma- 
trix Q can be parameterized as 

Q = Acxp(VF) , where AW + WA = (33) 

and A is the saddle point of the action (27) 



A(x,e) = 



1 2F(x,e) 
-1 



(34) 



The function F is related to the single particle distribu- 
tion function / through 

F(a:,e) = l-2/(a;,e). (35) 
The matrix W x e e r, in turn, is parameterized as follows: 



it" / Fx.e^x.e.e' ^x.e.e' H~ -^x.e^x.e.e' Fx.e' \ 

W x,e,e'-\ yj xt f , _«T, .TP . ) • ( 36 ) 



It is convenient to introduce the diffusion propagator 

+ D\7 2 )T>{x,x'u) = -l-5(x-x') . (37) 
ttTlv 
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The absence of diffusive motion in clean metallic leads 
implies that the diffusion propagator must vanish at the 
end points of the constriction. In addition, there is no 
current flowing in the transversal direction (hard wall 
boundary conditions). It follows that the component of 
the gradient of the diffusion propagator in that direction 
(calculated at the hard wall edges) must vanish as well. 
The correlation functions of the fields w, w are then given 
by: 



0.4 



0.3 



!^0.2 



(w(x,ei, e 2 )w(x', e 3 ,e 4 )) = 

2(27r) 2 <5(ei - e 4 )S(e 2 - e 3 )B{x, x', e l - e 2 ) , 

(w(x, ei,e 2 )w(x',e 3 , e 4 )) = -5(27r) 3 <5(ei - e 4 )<5(e 2 - e 3 ) 

(w(x, e 1 ,e 2 )w(x',e 3 , e 4 )) = . (38) 

To evaluate £3 one follows steps similar to those that 
led to the derivation of S2, see Ref. 17 . If all relevant en- 
ergy scales in the problem are smaller than the transver- 
sal Thouless energy (Eth = DjL\, where Lt is a width 
of a wire), the wire is effectively quasi-one dimensional. 
In that case only the lowest transversal mode of the dif- 
fusive propagator can be taken into account, which yields 



D(x 1 ,x 2 ) = — 



\Xl — X2 — X\ — X 2 + 



2xix 2 



(39) 



Here g ~ hvD. The electron distribution function in this 
system is equal to 



F(x,e) = -F eq [e 



eV 



eV 



(40) 



The quantities F and D determine the correlation func- 
tions, eq. (38). We can now begin to evaluate S3, (c.f. 
eq. (30)), performing a perturbative expansion in the 
fluctuations around the saddle point solution, eq. (34). 
After some algebra we find that in the zero frequency 
limit the third order cumulant is given by 



^ 3e 3 Airg 2 f L , , 
S 3 (uji = 0,uj 2 = 0) = — ^ — J dxidx 2 

/CO 
deF(e,a;i)D[0,a;i,i 2 ]V (F 2 (e,x 2 )) . (41) 
-00 



Integrating over energies and coordinates we obtain 



0.1 



-6-2 2 6 

log(eV/T) 

FIG. 2. The scaling function y plotted on a logarithmic 
scale, cf. eq. (42) 

Let us now discuss the main features of the function S3. 
In agreement with symmetry requirements S3 is an odd 
function of the voltage (the even correlator S2 is propor- 
tional to the absolute value of voltage), and vanishes at 
equilibrium. The zero temperature result (high voltage 
limit) has already been obtained by means of the scatter- 
ing states approach for single-channel systems 4 , and later 
generalized by means of Random Matrix Theory (RMT) 
to multi-channel systems (chaotic and diffusive) 18 . In our 
derivation we do not assume the applicability of RMT 19 . 
Our result covers the whole temperature range. We ob- 
tain that at low temperatures the third order cumulant 
is linear in the voltage 



(43) 



At high temperatures the electrons in the reservoirs are 
not anymore in the ground state, so the correlations are 
partially washed out by thermal fluctuations. One may 
then expand y(p), eq. (42), in a scries of eV/2T. The 
leading term in this high temperature expansion is linear 
in the voltage 



S 3 = -rl ■ 



(44) 



S 3 {u>i = 0,uj 2 = 0) = e 2 Iy(p) , 

= 6(-l + <*) + (! -26^ + ^ (42) 

where p = eV/2T. The function y is depicted in Fig.l, 
where it is plotted on a logarithmic scale. 



Note that although thermal fluctuations enhance the 
noise (compared with the zero temperature limit), 
eqs.(43) and (44) differ only by a numerical factor. The 
experimental study of S3 (and higher odd cumulants) 
provides one with a direct probe of non-equilibrium be- 
havior, not masked by equilibrium thermal fluctuations. 
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B. Strong Inelastic Collisions 

In our analysis so far we have completely ignored in- 
elastic collisions among the electrons. This procedure 
is well justified provided that the inelastic length greatly 
exceeds the system's size. However, if this is not the case, 
different analysis is called for. To understand why inelas- 
tic collisions do matter for current fluctuations, we would 
like to recall the analysis of 52 for a similar problem. The 
latter function is fully determined by the effective elec- 
tron temperature. Collisions among electrons, which are 
subject to an external bias, increase the temperature of 
those electrons. This, in turn, leads to the enhancement 
of 52, cf. Rcfs. 20,21 . In the limit of short inelastic length 

^c-c < L , (45) 

the zero frequency and zero temperature noise is 

5 2 (0) = ^-el . (46) 

In the present section we consider the effect of in- 
elastic electron collisions on 53. We assume that the 
electron-phonon collision length is large, l c - P h 3> L, 
hence electron-phonon scattering may be neglected. The 
Hamiltonian we are concerned with is 



H — H n + Hi nt . 



(47) 



The Coulomb interaction among the electrons is de- 
scribed by 



flint = ±J rfrrfr'*(r)*(r')^ (r-r')*(r)*(r') , 



where 



V (r-r') = 



r — r 



(48) 



(49) 



We need to deal with the effect of electron-electron inter- 
actions in the presence of disorder and away from equi- 
librium. Following Ref. 23 one may introduce an auxiliary 
bosonic field 



(50) 



which decouples the interaction in the particle-hole chan- 
nel. Now the partition function (eq. (25)) is a functional 
integral over both the bosonic fields Q and <j) , 



(Z) 



Q 2 



VQV(t>exp{iS total) ■ 



The action is 

iStotai = iS[$] + iS{<t>, Q] , 
= iTr{<P T VV$} , 

iS[$,Q] = -^Tr{Q 2 }+Trln 
4r 



(51) 



(52) 
(53) 

• (54) 



It is convenient to perform a "gauge transformation" 23 
to a new field Q 

Qt,t'(x) = exp(ik a (x,t)-y a )Qt,t'{x)exp(-ik a (x,t')j a ) . 

(55) 

Introducing the long derivative 

d x Q = VQ + i[Vk al a ,Q], (56) 
one may write the gradient expansion of eq. (54) as 



iS[Q, $] = w/Tr{($ - iujK) T l2 {<S> + i^K)} 



(57) 



7T1S 

T 



DTr{0 x Q} 2 + 4zTr{(e + {<f> a + iuk a )j a )Q} 



At this point the vector K T = (ki,k 2 ) that determines 
the transformation (55) is arbitrary. The saddle point 
equation for Q of the action (57) is given by the follow- 
ing equation 

Dd x {Qd x Q) + i[(e + {<j> a + iLuk a )~f a ), Q} = 0. (58) 
Let us now choose the parameterization 

Q = Aexp(I40, (59) 
where W represents fluctuation around the saddle-point 

1 2F[<j>](x,e)' 



A(x,e) = 







(60) 



Eq. (60) implies that the solution of the saddle point 
equation (58), determines F as a functional of <j>. We 
do not know, though, how to solve it. Instead we av- 
erage over (j) the eq.(58). The solution of this averaged 
equation, denoted by F, is determined by: 



(61) 



-DV 2 F{e) = I ee {F} , 
where the r.h.s. is given by 

I ee {F} = D J ^[{Vk 1 {u)Vk 1 {-Lo)){F{e)-F{e-Lo)) + 
{{Vk\u)Vk 2 {-Lo))-{Vk 2 {u)Vk\-Lo))){F{e)F{e^u)-l)]. (62) 

Since (A) is not a genuine saddle point of the action 
there is coupling between the fields W and <p, (Vfe) in 
the quadratic part of the action. However the coupling 
constant between those fields is proportional to the gradi- 
ent of the distribution (cf. eq.(84)). Therefore this term 
can be treated as a small perturbation. 

Taking variation of the action with respect to w, w, we 
obtain the following gauge, determining k[<j>]: 



DV 2 k 



DV 2 h +(f)i+ iuki = 2B[w, x] V 2 k 2 . 



2 — <P2 — i^k 2 = 



(63) 
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where 

B[w, x] = J de[l - F(e, x)F(e ~ uj, x)] . (64) 

Though we have failed to find the true saddle point the 
linear part of the action expanded around (A) is zero. It 
is remarkable to notice that under conditions (63) eq.(61) 
becomes a quantum kinetic equation 24 with the collision 
integral being I ee {F}. Coming back to our calculations 



we note that the correlation function of current fluctu- 
ations is a gauged invariant quantity (does not depend 
on the position of the Fermi level). This means that mo- 
menta q < \JujjD do not contribute to such a quantity 22 . 
In this case the Coulomb propagator is universal, i.e. 
does not depend on the electron charge. The fact that we 
address gauge invariant quantities allows us to represent 
the generating functional Z in terms of the fields Q and 
Vfc (rather than Q and <fi), as in Ref 23 . 



(Z)= JvVK cxp {-ivDTri^^V-^K}) J VQ exp ^ iS t [Q, VK] 



(65) 



Here we define 



D -i = l -DVi + iu6 x , x > 

! -£>V 2 — iu6 XiX > -2koS x<x iB u {x) 

where the expansion S = S° + S 1 + S 2 , is in powers of VK; the I — th power (I = 0, 1, 2) is given by 

iS°[Q] = ~ [dTt{VQ} 2 + 4i Tr{eQ} 



(66) 



(67) 



iS 1 [Q,VK] = -inu DTr{Vk al a QVQ} + Tv{{4> a +iu J k a ) 1 a Q} , 



(68) 



ttvD r 



Tr{Vfc a7 Q gVfc / 3 7 ' 3 Q}-Tr{Vfc Q 7 Q AVfc / 37' 3 A} 



From cq.(66) we obtain the gauge field correlation function 



(69) 



where 



Y(ui, x, x 1 ) 



(Vk a (x,u)Vki3(x' ,-u)) = —Y a ,p(u,x,x'), 



-2i-KVLo J dxi~D\-Ld, x, xi]B[u), xi]D[u), x\,x'] D[— w, x, x' 
D[u!,x,x'} 



Using eqs. (70,71) we rewrite eq.(62) for the quasi-one-dimensional wire as: 
DW 2 F(e) =I ee {e,x}, where 



I ee (e,x) = y J duj[ - 2i{jjTTvD[x,xi,-uj]B[uj,xi]D[xi,x,(jj]<yF{e) - F(e - w)) + 



(D[x,x,lu] -T>[x,x,-u])(l-F(e)F(e + uj))]. 



(70) 

(71) 

(72) 
(73) 



The total number of particles and the total energy of 
the electrons are both preserved during electron-electron 
and elastic electron-impurity scattering. The collision in- 
tegral, cq.(73), satisfies then 

/oc 
l ee (e,x)de = 0, (74) 
-OC 

/OC 
el ee (e,x)de = 0. (75) 
-oo 



We now consider the limit l ee L. The solution of 
eq.(61) assumes then the form of a quasi-cquilibrium 
single-particle distribution function 

fM)=lanh (iz_*r). (76) 

Here e is the total energy of the electron, and e<j> is 
the electrostatic potential and T{x) is the effective local 
temperature of the electron gas. To find the electrostatic 



8 



potential <f> we employ cq. (74). To facilitate our cal- 
culations we further assume that conductance band is 
symmetric about the Fermi energy and that the spectral 
density of single-electron energy levels is constant. Inte- 
gration over the energy, eq.(74) yields 



d x 2 (j>{x) = 0. 



(77) 



Solving eq.(77) under the condition that the voltage dif- 
ference at the edges of a constriction is V, we find 



ecj)(x) 



(78) 



Multiplying eq. (72) by energy and and employing eq.(75) 
we obtain an equation 



0. 



(79) 



The boundary condition of eq.(79) is determined by the 
temperature of the electrons in the reservoirs. Combin- 
ing eqs.((78) and (79)) we find the electron temperature 
in two opposite limits: 



T(x) 



^y/x(L-x) eV»T, 
eV^T . 



(80) 



Sa(ti — t 2 , t 2 — t-i) 



' -M(xi,ti)l (x 2 ,x 3 ,t 2 ,t 3 ) + (a;i,ti <-» 



(xuh ^x 2 ,t 2 ) + ^^-M(x 1 ,t 1 )M(X2,t 2 )M(xs,t 3 ) 



Eqs.((76), (78) and (80) determine the function F 
uniquely. We now replace the right-corner element of 
the matrix A (i.e. F[</>], cf. eq.(60)) by its average value 
F. 

To calculate S3 under conditions of strong electron- 
electron scattering (eq. (45)) one needs to replace the 

operators I and M in eq. (30) by their gauged values 

X 3 ,t 3 ) + 

(81) 



Vk,Q 



where the averaging is taken over the entire action S 
and the Gaussian weight function for Vif, as in eq. (51). 
Here we define (cf. eqs.(31), (32) with eqs. (82), (83)) 



I (x,x',t,t') = r mQ x ,t,t''y2Qx',t',t'Y2 -S t ,tHi \$x 



M(x, t) 



Tr^Jdh ([Qxjfyldx] Qx,t u 



72 



(82) 



(83) 



where the "long derivative", d x , is presented in eq.(56). 
In order to actually perform the functional integration 
over the matrix field Q we use the parameterization of 
eq. (59). We need to find the Gaussian fluctuations 
around the saddle point of the action (67,68,69). Though 
we did not find the exact saddle point, the expansion of 
Q around A works satisfactorily. The coupling between 
the fields Vfc and W which appears already in the Gaus- 
sian (quadratic) part is small, since it is proportional to 
the gradient of the distribution function: 



iS\ = -2i-ngTr^w x ,e,e> [Vfc 



lx,e' - 



e VF x . e - VF x ,,Vk 



lx,e' 



VF x ^Vk 2x , eU F x , e + F x ^Vk 2x . e ,. t VF x . e ] 



(84) 



Here the upper index refers to the power of the Vfc fields; 
the lower refers to the power oiw,w fields in the expan- 
sion. Considered as a small perturbation, iS\ does not 
affect the results. 

The more dramatic effect on the correlation function 
arises from the non-Gaussian part of the action, eqs. 



(68,69) (by this we mean non-Gaussian terms in either 
w,w or V-fQ. After integrating over the interaction an 
additional contribution to the Gaussian part (propor- 
tional to ww) of the action arises. To find the effective 
action zS^fW] we average over the interaction 16 . One 
notes that 



(85) 



(where, again, Sq refers to the component of the action, 
eq.(68), that has zero power of the w,w fields and one 
power of the Vfc field). In addition, due to the choice 
of the gauge, eq.(63), and the condition (l cc <C L), the 
averaging over Vfc does not generate terms linear in w, w 
in the effective action: 



-2inP I ^-<r,, A, /-' =0. (86) 



Combining eqs. (85 and 86) we find that the effective 
action acquires an additional contribution: 



cxp (1S 1 + iS 2 ) 



Vfc 



exp[{iS 2 ) + -{iS 2 iS 2 )) (87) 



The general form of the effective action is rather compli- 
cated, however for the low frequency noise only diagonal 
part of the action matters: 



TTV r 



iS 2 [w,w] = — Tr^w x 



DV 2 +l ee 



Wx.e.e ~ 



w x , e , e DVF; 



(88) 
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Here the operator 

i ee w x ,e,e = J duj[Yu(u))[w e , e -W e - U>e - U ] + 
(Y 12 (u)-Y 21 (w))[F e w e + F e - U w e>e ] + 



(89) 



is a linearized collision integral, i.e. a variation of the col- 
lision integral (73) with respect to the distribution func- 
tion. Substituting eqs.(82 and 83) into eq.(81) and cal- 
culating the Gaussian integrals with the action (88), we 
find 

^ 3e 3 Airg 2 f L , , 
63 (wi = 0,w 2 = 0) = — rT^5 — / dx\dx 2 



f 

J — ( 



de 1 de 2 F(e u x 1 )V[x u e 1 ;x 2 ,e 2 ]W (F 2 (e 2 ,x 2 )) , (90) 



where the "inelastic diffusion" propagator, T>, is the ker- 
nel the equation 

[-DW 2 - i ee ]V[ Xl , ei ; x 2 , e 2 ] = —6(6, - e 2 )S( Xl - x 2 ). (91) 

For weak electron-electron scattering the collision inte- 
gral is small, yielding the standard propagator of the 
diffusion equation, D (cf. eq. 39). In the presence of 
strong electron-electron scattering the collision integral 
dominates eq.(91). In this limit we evaluate the leading 
asymptotic behavior for (k n /L <C 1). On scales longer 
than the inelastic mean free path the distribution func- 
tion has a quasi-equilibrium form 

' e — n(x) — 5fi(x) " 



f(e,x) = f a 



T(x) + ST(x) 



where 



= ITS 



1 



(92) 



(93) 



+ cxp(a;) ' 

The values of the local temperature and electro-chemical 
potential can fluctuate. To find the correlations of these 
fluctuations we consider the equation (with the same Ker- 
nel as in eq.(91)) 

[-7JV 2 - 2 ee }6f(e, x) = SJ(e, x). (94) 

Integrating equation (94) with respect to energy and us- 
ing the particle-conservation property of the collision in- 
tegral (eqs.(74)) we find: 

H(x)= I dx'~D[x,x'] ( deJ(e,x'). (95) 

Using energy conservation (eq. 75) we find 

ST(x) = -^-^- I dx'B[x,x'] ( de(e- f i(x))SJ(e,x'). (96) 

7T T(X) J J 

Combining eqs. (95) and (96) we find: 

-nr 1 f 9 t ( ei ~M Xl ) 



D[xi,x 2 ](-1 



3 tx- ^,(xi)e 2 - ^(xx) 
tt 2 T( Xl ) T( Xl ) 



(97) 



Evaluating S3 explicitly we find that the third order 
current cumulant is 



5 3 (wi = 0,w 2 = 0) = 



D[xi,a;2] 



36e 3 Ag 2 eV 

T( Xl ) , 



> 3 Ag 2 eV f 1 

L 4 7T Jo 



dx\dx 2 



+ {xi - x 2 )- 



d 



T( Xl 



_T(x 2 ) ' T(x 2 ) dxi 

At high temperatures (cf. eq.80) one obtains 

5 3 (wi = 0,w 2 = 0) = Are 2 1, 
while at low temperatures 

5 3 (o; 1 =0, W2 = 0)=(^-^)e 2 J. 



(98) 



(99) 



(100) 



Our analysis was performed for a simple rectangular con- 
striction. However, our results hold for any shape of 
the constriction, provided it is quasi-one dimensional (we 
have considered a single transversal mode only). 



III. COMPARISON WITH THE KINETIC 
THEORY OF FLUCTUATIONS 

As was mentioned in Section I C the applications of 
the kinetic equation in the study of high order cumu- 
lants is not straight forward. To compare our microscopic 
quantum mechanical calculation with the semi-classical 
Boltzmann-Langevin scheme we consider the diagrams 
for the pair and third order correlation function, depicted 
mFig.(3). 

The diagram in Fig. (3-b) corresponds exactly to the 
results obtained above in the framework of the er-model 
formalism. Indeed, evaluating this diagrams, we recover 
the results (41) and (90). Comparing Fig. (3-b) with di- 
agram (3-a) (the latter determines the correlator (SJSJ) 
of Langcvin sources), we conclude that the third cumu- 
lant corresponding to Fig. (3-b) can be expressed in the 
form 



(57(1)57(2)57(3)) = 

,6(6 J (1)6 J (2)) 



dMA' 



W) 



7J>(4,4')(57(4')<5J(3)). (101) 



Indeed, the block on the right hand side of the diagram 
corresponds to the pair correlation function of random 
fluxes, while the left part can be obtained by functional 
differentiation of the diagram (3-a) with respect to the 
function F = 1 — 2 f. This justifies the regression scheme 
proposed recently by Nagaev 26 . Diagram (3-a) for the 
pair correlation function is local in space; the diagram 
for the third order correlator is not. Despite this non- 
locality the use of the Langevin equation (with non-local 
random flux) remains convenient; in that case the corre- 
lation function of the random fluxes needs to be calcu- 
lated from first principles. This is somewhat analogous 
to the description of mesoscopic fluctuations through the 
Langevin equation, proposed by Spivak and Zyuzin 25 . 
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(a) 



(b) 




FIG. 3. The second (a) and the third (b) current correlation functions. 
The vector vertices j represent the current operators, while collisions with static disorder correspond to the Hikami box (H). 
F = 1 — 2/, where / is the single-electron distribution function. D is the diffuson. 



Following Ref 17 the revisited version of the kinetic the- 
ory of fluctuations has been proposed by Nagaev 26 . He 
has noticed that for a diffusive system there exists a re- 
gression scheme of high order cumulants, expressed in 
terms of pair correlators 27 . However, a truly classical 
theory addressing high moments of noise needs to be 
taken on the same footing as the Boltzmann-Langevin 
theory, i.e. without appealing to quantum mechanical 
diagrammatics. Close inspection of the diagrams de- 
picted in Fig. 3 and the dimensionless parameter by 
which corrections are small implies that the regression 
recipe (eq. 101) has the same range of applicability as 
the Boltzmann-Langevin scheme itself. 

As a simple example of a classical problem for which 
the regression procedure can be applied, we consider the 
following simple model 6 . A heavy molecule of mass M 
and cross-section A is embedded in a gas of light clas- 
sical particles. The gas consists of N particles of mass 
to <C M, in thermal equilibrium with temperature T . It 
is enclosed in a narrow tube of volume V . The fluctuating 
velocity u of the molecule here is the analogue of the fluc- 
tuating electron distribution function Sf. The collisions 
of the molecule with the light particles is the counterpart 
of the electron's scattering on the random disorder. The 
motion of the molecule is governed by Newton law 



du , . 



(102) 



molecule. The latter can be calculated from 



F{t) = 2mA 



o 



dvv 2 P 1 {v + u{t),t)- dvv 2 P 2 (v+u(t),t) 



(103) 



where Pi and P 2 are the fluctuating distribution func- 
tions of particles on the left and on the right sides of the 
system respectively. Since the time between two consec- 
utive collisions is much shorter than the relaxation time 
of the molecule (r m ) one can average over fast collisions 
while considering the relaxation dynamics of the heavy 
molecule: 



du 
~db 



1 



-u. 



(104) 



To find the various velocity correlation functions, one 
needs to know the corresponding correlation functions of 
the random forces. Using the fact that the equilibrium 
fluctuations of the distribution function in a Boltzmann 
gas are Poissonian ( 7 , Section 114) and neglecting the 
small velocity u relatively to v, one finds the pair corre- 
lation function 



(SF(t)SF(t')) = 



16 AN 
V 



m 2 (v 2 ) 3/2 S(t-t'). (105) 



Eq. (105) yields the correct value for the thermal fluctu- 



Here u is the velocity and F is the force acting on the ations 



T/M. To find the higher order correlations 
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of the random force one has to take into account the de- 
pendence of the force on the velocity of the molecule. We 
expand eq.(103) up to second order in u. After averag- 
ing over all possible pair correlations (such as (uSP) and 
{5P5P}) we obtain 

(SFWSFfotfFfo)) = 



512 {AN\ 2 i, 



2tt V V J 



i»V)M« s(h-t 2 )9{t 2 -t 3 ) y (106) 



}■ 



Here p denotes all permutation over (ti,t2,ts). The 
triple force correlator was reduced to pair correlators 
only. Bearing in mind the analogy between u and Sf 
on one hand, and the scattering of the electrons on the 
disorder and the scattering of the molecule by light par- 
ticles on the other hand, we note that the reduction of 
the third order cumulant of the random forces to the pair 
correlation functions is similar in both cases. 

We finally come back to the question of whether it is 
possible to calculate high order cumulants employing the 
classical kinetic equation (Boltzmann-Langevin) rather 
than resorting to the diagrammatic reduction scheme de- 
picted above. To be able to answer this question we first 
note that the applicability of the kinetic theory requires 
that both l/g,h/T COII kBT <C 1. Here T cor is the correla- 
tion time for the current signal and g is the dimensionlcss 
conductance. For non-interacting electrons Onsager re- 
lation and Drude formula yield t coi — r, where r is the 
transport time. It follows that the condition 



n 



z = 



gT cor k B T 



< 1 



(107) 



must be satisfied. The reason for theses inequalities are 
first, that g needs to be large in order for quantum in- 
terference effects to be negligible. Secondly, the value 
h/k^Tr needs also to be small; the transport equation 
is applicable for times longer than the duration time of 
an individual collision St, (t >• St). For a degener- 
ate electron gas the uncertainty relation requires that 
St > h/k B T 2 ^ 29 . 

We now evaluate the relative magnitude of, e.g., the 
fourth and the second cumulants. We consider the ratio 



«/ 4 (0) » 
«/ 2 (0) » 2 



(108) 



where < < > > denotes the irreducible part of the corre- 
lator. At equilibrium 

z~z. (109) 

As we see, for the kinetic theory to be valid, the value of 
the parameter z needs to be small. But this is exactly 
the parameter (z) by which the high order correlation 
functions are smaller than the lower ones (cf. eq.107). In 
other words, the evaluation of high order cumulants goes 
beyond the validity of the standard Boltzmann-Langevin 



equation. This is why we have to resort to the diagram- 
matic approach: either to justify the reduction of high 
order cumulants to pair correlators, or, alternatively, to 
introduce non-local noise correlators in the Boltzmann- 
Langevin equation. 



IV. COUNTING STATISTICS 

So far we have studied the second and third order cu- 
mulants. In the present section we discuss the whole 
distribution function of the low frequency electron cur- 
rent (so-called counting statistics). The zero tempera- 
ture limit had been studied by Levitov et. al. 18 . The full 
temperature regime was addressed by Nazarov 31 . Here 
we present a different derivation based on the a-model 
approach. 

Being a stochastic process, the charge transmission can 
be characterized by the probability distribution function 
(PDF) Pt(n) of the probability for n electrons to pass 
through the constriction within the time window t. In 
practice it is more convenient to work with the Fourier 
transform of the PDF, the characteristic function 



(110) 



Below we evaluate kj[A] for the case of elastic scatter- 
ing. By expanding the logarithm of characteristic func- 
tion over its argument we can find the cumulants (irre- 
ducible correlation functions) of a transmitted charge 



k 



(111) 



For the problem of diffusive junction the disorder average 
characteristic function can be represented as: 

«f[A] = JvQexp(iS[Q,a}), (112) 
where the external source is given by 

*® = { f others . ( 113 ) 

The applied bias enters the problem through boundary 
conditions on Q at the edges of the constriction: 



A(0,e) 



,A(L,e) 



2F 

-1 



and F is defined by eq.(35). 

Inasmuch as we are not interested in spatial correla- 
tions, the external source term is a function of time only. 
Therefore, by performing the transformation 

Q{x,t,t') = e- ixa2(t) ^Q{x,t,t')e ixa2{t '^\ (114) 
one gets: 
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iS[Q,a\ = Tr\ d(vQ) - 



4r 



d_ d_ 

dt + dt' 



7 )Q(t,t') + 



i fJ^a 2 (t)^2Q(t,t') - Q(t,t')i-^a 2 (t'h 2 } , (115) 

and the boundary conditions change correspondingly: 
Q{0,t,t') = A(0,t-t'), 

Q{L,t,t') = e- W2(t)72 A(L,t-0e W2(t ' )72 , (116) 
with 

A(t-i') = / ^e 4e ^ t ')A(e). (117) 

J 27T 

We use the saddle point approximation to calculate the 
characteristic function, Eq. (112). In the presence of an 
external potential the minimum of the action, cq. (115), 
satisfies 

dv ( qvq)-(| + | 7 )q + 

d d 
i-g- t a 2 (th 2 Q(t, t') - Q(t, t')i—a 2 {t') l2 = 0. (118) 

Let us define a parameter that can roughly be re- 
garded as "the number of attempts per channel" , M — 
maxjeV, T}t/h. We will focus on the case where the 
time window is much larger than the Thouless time 
i S> txhouiess and the conductance g 3> M 3> 1. In- 
side the region < t, t' < t eq. (118) becomes (up to the 
corrections 0(t T hou\css/t)) 



DV{QVQ) = 0. 



(119) 



This can be represented in the form resembling a current 
conservation law 



DVJ = 0, 

where the current is defined as 



(120) 



J=QVQ. (121) 

The solution of Eq. (120), Q sp (x) can be written as 

g sp = g(0)cxp(Ja;), (122) 

and from the boundary conditions (eq. 116) it follows 
that 



J = HQ(0)Q(L)). 



(123) 



One may show that the anticomutator of the matrix Q 
with J vanishes 



{g(o),j} + = o. 



(124) 



Based on this property, we can show that the following 
statements do hold: 

(1) The square of the matrix Q is still equal to unity 

Q%{x) = 1, (125) 

(2) The determinant of the matrix Q satisfies 

Det[g sp (x)] = -1, (126) 

(3) The trace of the matrix Q is equal to zero 

Tr{Q sp (x)} = 0, (127) 

(4) The trace of the matrix J is zero 

Tr{ J} = 0, (128) 

(5) The square of the matrix J is proportional to the unit 
matrix 



J 2 = S[e,X]I, 



(129) 



where 5[e, A] is given by (there is small discrepancy with 
the result obtained in Ref. 31 ) 



<5[e,A] =ln 



(2/i - l)(2/ 2 - 1) + 2.M1 - h)^ + 2/ 2 (l - fi)e + 



2 v /(e iA/2 - e-^)(l - h + he-^)(l - f 2 + f 2 e^)(f 2 (l - h)e^ - h (l - / 2 ) e -<V2) 



(130) 



Using these properties we find the disorder averaged 
counting statistics 



k[\] = exp 



Gt 



deS 2 [e,X\ . 



(131) 



S 2 
S 4 



2GT 

- e 2 GT. 
3 



(132) 



The zero temperature limit coincides with the results de- It is worthwhile to note that the value of £4 agrees with 
rived previously by Lee et. al. 18 . the one obtained in Ref. 26 . As we see the counting statis- 

At finite temperatures we find: tics of the current in a disordered wire is not Gaussian. 
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Remarkably, all the information contained in the count- 
ing statistics can be extracted from the pair correlation 
function (of the distribution function), eq.(17). 

Finally we would like to discuss the role of inelas- 
tic electron-phonon scattering. As has been already 
realized 1 , such an interaction suppresses shot noise. 
Moreover, based on our approach, one can show that 
in the limit of ^ _ p h <C L (macroscopic conductor) the 
current fluctuations are Gaussian (to leading order in 
L/le-ph)- To show this, we repeat our analysis concern- 
ing S3. We find that eq. (90) still holds, but the dis- 
tribution function and the inelastic diffusion propagator 
need to be calculated taking electron-phonon collisions 
into account. In the presence of both electron-electron 
and electron-phonon interaction one obtains 



[dV 2 f + i ee + r- ph ]F = 



(133) 



(cf. eq.61) where X e ~ ph is the linearized electron-phonon 
collision integral. The inelastic diffuson (cf. eq.91) is now 
determined by 

hDV 2 -f ee -i e -P h ]V[x l7 e 1 ;x 2 ,e 2 ] = 
—5(e 1 -e 2 )5(x 1 -x 2 ). (134) 

In the limit L/7 c _ ph 3> 1 fluctuations of the chemical 
potential is the only long-range propagating mode in 
the problem (no fluctuations of the fcgT along the sys- 
tem). Solving eqs.(133, 134 and 90) we find that S3 
vanishes as (l e ^ p h/L) 2 . A conductor longer than the 
electron-phonon length (Z c -ph) can be viewed as a num- 
ber (L/Z c _ph) of resistors connected in series. The cur- 
rent fluctuations of a single resistor are given by eq.(131) 
and would render the corresponding voltage fluctuations. 
Since the fluctuations at different resistors are uncorre- 
cted (local fluctuations in the chemical potential do not 
change the resistance significantly) the large number of 
the resistors results in Gaussian current fluctuations. By 
contrast, temperature fluctuations (in the case of solely 
electron-electron interactions) have long distance corre- 
lations, and give rise to non-Gaussian fluctuations. 
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